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分 数 阶 模型 模拟 混凝土 徐 变 特性 的 多 参数 识别 
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摘 要 :本 研究 主要 研究 混凝土 徐 变 特性 的 分 数 阶 模型 拟 合 及 其 多 参数 识别 问题 。 在 模型 方面 , 提 
出 利用 修正 的 分 数 阶 Maxwell 模型 和 分 数 阶 Poynting-Thomson 模型 两 种 模型 模拟 混凝土 的 徐 变 实 
验 数 据 , 并 将 结果 进行 对 比 ,通过 数据 拟 合 和 误差 分 析 验 证 两 种 模型 的 有 效 性 。 在 多 参数 识别 方 
面 ,分 别 采 用 贝 叶 斯 算法 和 布谷 鸟 搜索 算法 两 种 算法 识别 两 个 模型 中 的 多 个 参数 。 研 究 表 明 : 修 正 
的 分 数 阶 Maxwell 模型 和 分 数 阶 Poynting-Thomson 模型 在 刻画 混凝土 徐 变 特 性 中 均 是 有 效 的 ; 贝 叶 
斯 算法 和 布谷 鸟 搜索 算法 在 分 数 阶 模型 的 多 参数 估计 问题 中 均 是 可 行 的 ,但 布谷 岛 算 法 搜索 速度 
更 快 . 误 差 更 小 效率 更 高 ,在 分 数 阶 模型 的 多 参数 识别 问题 中 性 能 更 优 。 
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Fractional-order model simulation and multi-parameter 
identification for the creep behavior of concrete 
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(School of Science, Jiangnan University ,214122 Wuxi ,China ) 


Abstract :Focusing on the creep behavior of concrete ,the fractional-order model simulation and the corre- 
sponding multi-parameter identification are studied and analyzed in this paper. Both the modified fractional 
Maxwell model and the fractional Poynting-Thomson model are proposed to simulate experimental data of 
concrete creep. The results are compared ,and the validity of the two models is verified by data fitting and 
error analysis. As for the corresponding multi-parameter identification ,both the Bayesian algorithm and the 
cuckoo search algorithm are used to identify the multiple parameters in the two models. Results show that 
the modified fractional Maxwell model and fractional Poynting-Thomson model are both effective in descri- 
bing the creep characteristics of concrete. Both the Bayesian algorithm and the cuckoo search algorithm are 
feasible in the fractional multi-parameter estimation problem ,but the cuckoo algorithm has higher search 
speed , smaller error ,higher efficiency ,and better performance in the multi-parameter identification problem 
of the fractional order model. 
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混凝土 是 当代 最 主要 的 建筑 材料 之 一 ,被 广泛 
应 用 于 土木 工程 .造船 业 .地 热 工 程 等 行业 。 研 究 混 
凝 土 的 徐 变 ,松弛 等 特性 ,对 现代 社会 的 建筑 稳定 性 
和 安全 性 具有 重要 意义 。 
近年 来 ,已 有 不 少 学 者 致力 于 研究 混凝土 徐 变 
特性 的 本 构 模 型 表达 。 由 于 整数 阶 导数 算 子 是 
局 部 算 子 ,不 能 很 好 地 描述 非 局 部 行为 ,而 分 数 阶 导 
数 算 子 是 非 局 部 算 子 ,具有 时 间 上 的 记忆 效应 和 空 
间 上 的 长 程 依赖 性 。 近 年 来 ,分 数 阶 本 构 模 型 成 为 
秋 弹 性 材料 研究 领域 的 热点 之 一 ”。 张 为 民 等 中 
研究 表明 ,利用 分 数 阶 导数 模型 模拟 混凝土 的 黏 弹 
性 特性 ,误差 更 小 ,计算 效果 更 好 。 梅 生 启 “构建 了 
一 种 修正 的 分 数 阶 Maxwell 模型 ,模拟 了 混凝土 的 
徐 变 及 其 徐 变 恢复 特性 ,并 利用 非 线 性 最 小 二 乘 算 
法 给 出 了 模型 中 的 参数 取 值 。 朱 端 等 “ 基于 分 数 阶 
Burgers 模型 研究 了 早 龄 期 混凝土 的 徐 变 特性 。 
与 分 数 阶 导 数 建 模 的 实际 应 用 相关 联 的 一 个 重 
要 问题 是 确定 模型 中 的 重要 参数 "” 。FAN 等 (中 利 
用 贝 叶 斯 算法 研究 了 广义 分 数 阶 单元 网 络 Zener 模 
型 的 参数 估计 问题 。 杨 秀 '"! 对比 分析 了 拟 牛 顿 算 
法 . 非 线性 共 斩 梯 度 方 法 与 贝 叶 斯 算法 在 分 数 阶 模 
型 参数 估计 问题 中 的 可 行 性 ,结果 表明 ,对 于 多 参数 
识别 问题 , 贝 叶 斯 算法 更 加 有 效 。CHI 等 利用 布 
谷 乌 算法 实现 了 对 分 数 阶 微分 方程 的 参数 估计 。 
鉴于 目前 对 混凝土 材料 的 研究 大 多 集中 于 模型 
拟 合 ,而 对 相关 模型 的 多 参数 识别 问题 研究 较 少 ,本 
研究 基于 混凝土 材料 的 徐 变 实验 数据 ,研究 分 数 阶 
模型 在 模拟 混凝土 黏 弹 特 性 中 的 应 用 及 其 多 参数 识 
别 问题 。 提 出 利用 修正 的 分 数 阶 Maxwell 模型 和 分 
数 阶 Poynting-Thomson 模型 模拟 混凝土 的 徐 变 数 
据 , 并 将 研究 结果 进行 对 比 , 分 析 两 种 分 数 阶 模型 在 
模拟 混凝土 徐 变 特性 中 的 可 行 性 ;针对 分 数 阶 本 构 
模型 ,研究 相关 的 多 参数 识别 问题 ,分 别 利 用 贝 叶 斯 
算法 和 布谷 鸟 搜 索 算 法 同时 识别 分 数 阶 模型 中 的 多 
个 参数 。 最 后 给 出 数值 算 例 ,分 析 贝 叶 斯 算法 和 布 
谷 鸟 搜 索 算 法 在 混凝土 分 数 阶 徐 变 模型 的 多 参数 识 
别 中 的 有 效 性 。 


1 分 数 阶 锋 弹性 本 构 模 型 
本 研究 主要 研究 修正 的 分 数 阶 Maxwell 模型 和 


分 数 阶 Poynting-Thomson ( FPT) 模型 在 模拟 混凝土 
徐 变 特性 方面 的 可 行 性 。 
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1.1 修正 的 分 数 阶 Maxwell 模型 


分 数 阶 Maxwell 模型 由 一 个 弹性 元 件 及 一 个 分 
数 阶 元 件 串联 而 成 ,而 修正 的 分 数 阶 Maxwell 模型 
则 是 在 分 数 阶 Maxwell 模型 的 基础 上 增加 一 个 黏 性 
元 件 , 如 图 1 所 示 。 


图 1 修正 的 分 数 阶 Maxwell 模型 
Bb Wonied naotional Maewell od 
修正 的 分 数 阶 Maxwell 模型 的 本 构 关 系 为 
l-a 

> | We ; TI De a 9 
其 中 : go (i) 为 应 力图 数 ; e(b 为 应 变 函 数 ;表示 时 
间 ; 五 为 弹性 系数 ; Ai 为 分 数 阶 元 件 的 单元 参数 ; ri 
为 松弛 时 间 ; 0 < a 1 为 分 数 阶 导数 阶 数 ;jw 为 竹 
性 系数 ; 4 所 DZXdz 为 连续 孔 数 fi1) 的 Riemann- 
Liouville 导 数 , 且 


df(t 1 dr 了 
全 FTFG-a) qd 上 
(tt >0,0 <a<1) 
其 中 , 研 (…) 为 Gamma 因数 。 
在 荷载 持续 作用 下 ,修正 的 分 数 阶 Maxwell 模 
型 ”的 徐 变 柔 量 可 得 为 
1 t t 
A 中 
当 持 傈 oo 时 ,修正 的 分 数 阶 Maxwell 模型 在 持 
谷 阶段 的 应 变 表 达 式 为 
a(t) = J(t) "oo+ eo 
1 1° 


= | + - 人 
E prmi*:T(l+ta) kW 


CQ 


) au + Eo 
(3) 
其 中 , su 为 初始 应 变 。 
1.2 分 数 阶 Poynting-Thomson 模型 
FPT 模型 由 一 个 弹性 元 件 .一 个 黏 性 元 件 和 一 
个 分 数 阶 元 件 “构成 ,如 图 2 所 示 。 


FPT 模型 的 应 力 应 变 本 构 关系 为 
Wa do(D) MT, dog() 


EO FE dt FE dz 
de(t) a d"e(t) 
一 HK» di 十 Wi171 d 尖 (4) 
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图 2 FPT 模型 
Fig.2 上 PT model 
在 荷载 持续 作用 下 ,FPT 模型 的 徐 变 柔 量 悟 ] 为 
1 上 在 -7 
J(i) dl + (5) 
i Wa/ M1 
其 中 已 ， 为 广义 Mittag-Leffler( M-L) 函数 ,定义 为 


Z 
Ez2)= 2 Fm 


,(zEeEC,p>0,g >0) 
f=0 + g) 2 . 


当 持 荷 为 oo 时 ,FPT 模型 在 持 荷 阶段 的 应 变 表 
e(t) = J(1) :oo + eo 


l-a a 
ee 7i 
Ek wk | M2/ WI 


:mm 十 CE0 


(6) 
当 FPT 模型 中 的 弹性 元 件 和 黏 性 元 件 蔡 换 为 分 
数 阶 单元 时 ,模型 推广 为 广义 的 分 数 阶 FPT 模 
型 。JIANG 等 “研究 表明 分 数 阶 FPT 模型 比 一 
般 的 组 合 模型 拟 合 精度 更 高 , 既 能 避免 整数 阶 模型 
在 里 变 曲线 扬 点 附近 与 试验 数据 拟 合 效果 差 的 浆 
端 ,又 能 反映 蠕 变 的 非 线性 渐变 过 程 ,可 以 更 好 地 描 
述 材料 的 竹 弹 特性 。 


2 多 参数 识别 算法 


2.1 贝 叶 斯 算法 


贝 叶 斯 算法 的 基本 原理 在 于 利用 与 未 知 变量 有 
关 的 统计 数据 来 获取 未 知 变量 的 条 件 分 布 函数 ” 。 
在 先 验 信息 较 少 的 情况 下 , 相 较 于 极 大 似 然 估计 法 ， 
贝 叶 斯 算法 的 模拟 效果 更 佳 ”。 贝 叶 斯 基本 公 
式 为 


i _ p(y| 0)p(0) (7) 
p(y) 


其 中 : 9 代表 待 估计 参数 ; y 代表 与 之 相关 的 统计 数 
据 ; p(01y) 代表 后 验 概率 密度 函数 ; p(y1 0) 为 似 
然 函 数 ; p(9) 是 未 知 变量 9 的 先 验 概 率 密 度 函 数 。 
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在 实验 数据 y 是 已 知 量 , 则 贝 叶 斯 定理 也 可 以 表 
示 ? 


PCOI7Y) «p(yl| 0)p(0) (8) 
基于 后 验 概率 密度 函数 ,可 采用 多 种 方法 进行 


未 知 量 的 点 估计 ,如 最 大 后 验 估计 (MAP): 6, = 


aug maxop( 91 y) ,或 后 验 均 值 估计 :x = Ey,01y。 

本 研究 采用 后 验 均值 估计 对 模型 中 的 未 知 参 数 
进行 识别 。 对 于 待 识别 参数 的 先 验 概率 分 布 ,本 研 
究 均 选 为 高 斯 分 布 。 假 设 参数 9 服从 均值 为 u , 方 
差 为 v 的 高 斯 分 布 , 记 作 9 ~ N(u,v,) , 则 6 的 先 验 
概率 分 布 函数 为 


/(0) = 


oO 


此 时 , 似 然 函数 p(y1 0) 表达 式 为 
p(y1 60) = 


yexp( - (se(0) 0 3 (10) 
(270, ) 3 Zu 
其 中 : s(0) 为 参数 为 6 时 的 模型 计算 值 ; v。 为 统计 
数据 y 中 独立 同 分 布 的 高 斯 随机 噪声 的 方差 ; 入 为 
统计 数据 y 的 长 度 。 
贝 叶 斯 算法 中 ,从 当前 状态 过 渡 至 新 的 状态 时 
设 定 了 一 个 接受 概率 ,具体 表示 为 
Pb) .p(y1 0;) .gq(0;1 0.) 
"p(0;) POLO) g(0,1 0.) 
(11) 
其 中 : 9, 表示 当前 状态 的 第 i 个 参数 的 值 ; 9; 表示 建 
议 的 新 状态 的 第 i 个 参数 的 值 ; g( 0; 1 0;) 为 建议 分 
布 函 数 , 一 般 表达 式 为 
q(0;| 0;) ~ N(0;,0,) (12) 
贝 叶 斯 算法 的 算法 流程 如 算法 1 所 示 。 
算法 1 贝 叶 斯 算法 
1) 确定 参数 的 先 验 分 布 p(01),p(9,),…， 
pP(0,) 及 参数 的 初始 值 ; 
2)whilej < 最 大 允许 迭代 次 数 ，; 
3)foris =1:n; 
4) 产 生 服从 建议 分 布 g( 9,1 9,) 的 新 参数 0 ; 
5 ) 计算 接受 新 参数 的 接受 概率 P, ; 
6) 产 生 [0,1] 之 间 的 随机 数 > 与 P, 比较 判断 ， 
若 r < P,; 则 接受 新 参数 0; ,否则 保持 原 值 ; 
7 ) end for; 
8 ) end while。 


P, = min{l 
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2.2 布谷 鸟 搜索 算法 


布谷 乌 搜 索 算 法 (cuckoo search, CS ) 是 一 种 群 
智能 优化 算法 ,常用 于 求 优化 问题 的 近似 解 。 它 模 
拟 布谷 乌 的 巢 寄 生殖 行为 寻找 最 优 解 ,同时 利用 莱 
维 飞行 确定 布谷 鸟 每 次 搜索 的 步 长 ”。 为 了 便于 
描述 布谷 乌 的 梨 寄 生殖 行为 ,在 布谷 鸟 搜 索 算法 中 
有 以 下 3 个 理想 规则 。 规 则 :每 只 布谷 鸟 一 次 只 
产 一 个 卵 ,并 随机 选择 一 个 位 置 的 鸟 代 进行 孵化 ;在 
随机 选择 的 一 组 鸟巢 中 , 质 优 的 马 巢 将 会 被 保留 到 
下 一 代 ; 可 利用 的 马 巢 数量 固定 ,鸟巢 主人 发 现 外 来 
鸟 蛋 的 概率 为 Pe [10,1] 。 当 鸟 集 主人 发 现 外 来 的 
布谷 岛 蛋 的 时 候 , 它 会 将 布谷 岛 蛋 丢弃 或 者 重新 建 


立新 的 鸟巢 。 
基于 上 述 规 则 ,布谷 鸟 位 置 更 新 公式 为 
X= +yL(s,A) (13) 


其 中 : XX; 表示 第 i 个 乌 划 在 第 1 代 的 鸟巢 位 置 ; @ 表 
示 点 对 点 乘法 ; y > 0 表示 步 长 的 控制 量 , 取 为 向 
量 ,具体 计算 为 y = yo( 忆 一下) ,yo 为 常数 ,XX 、 总 
为 种 群 中 随机 选择 的 两 个 个 体 ; 2(s,A) 为 莱 维 飞行 
的 随机 搜索 路 径 ,其 莱 维 飞行 的 方向 服从 均匀 分 布 ， 
其 行走 步 长 服从 莱 维 分 布 ” , 即 

ATCA)sin( 他) 


1 
(sy VE ym 


(s 0) 


(14) 
其 中 ; 1 < B < 2 ;J 和 ， 均 为 服从 正 态 分 布 的 随机 
数 ,pp ~ N(0,02) ,0 ~ N(0,0?) ,0, = 1 ,参数 wm 
的 计算 公式 为 ” 


T(1 +B)sin(mB/2) 1 
民 和 全， | (> 

布谷 鸟 搜索 算法 的 算法 流程 如 下 。 
算法 2 布谷 乌 搜 索 算法 

1 ) 初 始 化 种 群 ,设置 P=0.25，; 

2)while g < 最 大 迭代 次 数 ; 

3 ) 计 算 每 个 岛 巢 的 适应 度 ; 

4) 保 留 最 优 的 乌 站 位置, 每 次 迭代 都 对 其 进行 
更 新 ; 

5 ) 利 用 莱 维 飞行 得 到 新 的 位 置 ; 

6) 与 上 一 代 乌 巢 比 较 , 择 优 保存 ; 

7) 产 生 [0,1] 之 间 的 随机 数 r 与 P 比较 判断 ， 
若 r < P 则 接受 新 鸟 丰 ,否则 保持 原 值 ; 

8 ) end while。 
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3 混凝土 徐 变 实验 数据 拟 合 


为 了 验证 修正 的 分 数 阶 Maxwell 模型 和 FPT 模 
型 在 模拟 混凝土 徐 变 特性 中 的 可 行 性 及 参数 识别 算 
法 的 有 效 性 ,本 人 研究 分 别 采 用 修正 的 分 数 阶 Maxwell 
模型 和 FPT 模型 两 种 模型 对 C40 强度 混凝土 在 荷 
载 速率 为 1 kN/s 时 的 徐 变 测试 数据 ”进行 拟 合 。 
在 参数 识别 过 程 中 ,分 别 采 用 贝 叶 斯 算法 和 布谷 乌 
算法 确定 分 数 阶 模型 中 的 多 个 参数 ,并 将 两 种 算法 
所 得 结果 进行 对 比 ,分 析 两 种 算法 的 有 效 性 。 

为 了 刻画 分 数 阶 模型 对 实验 数据 的 拟 合 效果 ， 
记 模 型 模拟 结果 与 实验 数据 之 间 的 误差 为 


N 
6 洁 3 [log(e™™ (ti) +1) -log(e™ (ti) +1)1 
(16) 
其 中 : e™(1) 表示 模型 计算 值 ; ae (1) 表示 实验 数 
据 测 量 值 。 在 布谷 马 算 法 中 ,采用 模型 输出 与 实验 
数据 的 误差 倒数 作为 适应 度 对 算法 结果 进行 评价 ， 
即 适应 度 为 1/e 。 

在 数值 算 例 中 ,弹性 系数 瑟 取 为 32.39 GPa。 在 
修正 的 分 数 阶 Maxwell 模型 中 ,得 识别 参数 为 9，= 
(a nprn) ,在 FPT 模型 中 得 识别 参数 为 
0 = (QW AT ) o 

对 于 修正 的 分 数 阶 Maxwell 模型 和 FPT 模型 ， 
模型 中 的 分 数 阶 单元 由 3 个 参数 ( w,w ,Ti ) 表征 ， 
然而 ,通过 数值 实验 发 现 ,虽然 参数 jw，、7| 的 取 值 不 
稳定 ,但 jw7? 始终 趋 于 一 个 稳定 值 , 故 在 参数 估计 
结果 中 ,将 7? 作为 一 个 整体 来 分 析 。 


3.1 修正 的 分 数 阶 Maxwell 模型 拟 合 及 参数 识别 


利用 修正 的 分 数 阶 Maxwell 模型 拟 合 混凝土 的 
徐 变 实验 数据 ,并 确定 模型 参数 。 在 贝 叶 斯 算法 中 ， 
参数 初始 迭代 值 0 = (op rm) 取 为 (0.5， 
1,1,1) ,在 布谷 岛 算法 中 ,参数 初 值 在 ay e [0,1] ， 
Am e[0,20] ,wo es[0,5] , 71 e [0,100] 范围 内 
随机 生成 ,参数 B 取 为 1.5。 

图 3 为 利用 贝 叶 斯 算法 所 得 的 修正 的 分 数 阶 
Maxwell 模型 中 参数 gat、1/(ww7T%) 、1Ms 的 迭代 
结果 图 。 图 4 为 利用 布谷 鸟 搜索 算法 所 得 的 参数 
ma JIZ(arn) 1 的 迭代 结果 图 。 由 图 可 见 ,在 
迭代 一 定 次 数 后 ,两 种 算法 所 得 的 参数 均 趋 于 稳定 
值 , 说 明 贝 叶 斯 算法 和 布谷 鸟 搜 索 算 法 在 分 数 阶 模 
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型 参数 估计 问题 中 是 均 是 收敛 的 。 但 布谷 乌 搜 索 算 
法 达到 稳定 时 所 需 的 迭代 次 数 更 少 , 说 明 布 谷 鸟 搜 
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ee 问题 中 效率 更 高 , 速 
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图 3 修正 的 分 数 阶 Maxwell 模型 贝 叶 斯 算法 迭代 结果 


Fig.3 The Bayesian algorithm iteration results of modified fractional Maxwell model 
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图 4 修正 的 分 数 阶 Maxwel 模型 布谷 鸟 搜索 算法 迭代 结 


Fig.4 The cuckoo search algorithm iteration results of modified fractional Maxwell model 


图 5 表示 布谷 鸟 搜索 算法 的 适应 度 迭 代 图 。 可 
以 看 出 , 随 着 迭代 次 数 的 增加 ,模型 输出 与 实验 数据 
的 误差 逐渐 减 小 且 趋 于 一 个 稳定 值 ,进一步 说 明 布 
谷 鸟 搜 索 算法 在 分 数 阶 多 参数 识别 中 是 有 效 的 。 
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办 代 次 数 
5 修正 的 分 数 阶 Maxwell 模型 布谷 鸟 算 法 适应 度 迭代 图 


Fig.3 The fitness iteration results of cuckoo search ns 


for modified fractional Maxwell model 
图 6 显示 了 修正 的 分 数 阶 Maxwell 模型 对 混 凝 
土 的 徐 变 实验 数据 的 拟 合 效果 。 其 中 , 实 线 和 虚线 
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分 别 表示 将 贝 叶 斯 算法 和 布谷 马 算法 所 得 的 参数 取 
值 带 入 模型 得 到 的 拟 合 效 果 。 由 图 6 可 见 , 两 条 曲 
线 与 实验 数据 完美 的 拟 合 在 一 起 ,说 明 贝 叶 斯 算法 
与 布谷 岛 算 法 所 得 的 参数 取 值 均 是 有 效 的 。 考 虑 到 
布谷 乌 算 法 的 搜索 速度 更 快 , 氨 代 次 数 更 少 , 布 谷 乌 
搜索 算法 在 该 分 数 阶 模型 的 多 参数 识别 问题 中 表现 
更 做 。 


口 “数据 
一 一 一 布谷 鸟 算法 
贝 叶 斯 算法 
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持续 时 间 /s 
6 修正 的 分 数 阶 Maxwell 模型 拟 合 实验 数据 
Fig.6 Experimental data fitting of modified 


fractional Maxwell model 
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表 1 为 修正 的 分 数 阶 Maxwell 模型 的 参数 佑 计 
结 
表 1 修正 的 分 数 阶 Maxwell 参数 识别 结果 
Tab.1 The parameter estimation results of modified 


fractional Maxwell model 


参数 贝 叶 斯 算法 布谷 鸟 算法 
a 0.61 0.63 
1X(Cwirg ) 0.040 5 0.038 1 
1/ 0. 000 363 0.000 018 3 
e 0. 000 385 0. 000 378 


可 以 看 出 ,两 种 算法 所 得 的 参数 估计 结果 稍 有 


偏差 ,但 分 析 误 差 。 可 见 ,虽然 参数 取 值 有 偏差 ,但 
均 与 实验 数据 拟 合 效 果 较 好 ,说 明 在 利用 修正 的 分 
数 阶 Maxwell 模型 拟 合 混凝土 的 徐 变 实验 数据 时 ， 
相应 的 多 参数 识别 问题 存在 解 不 唯一 的 情形 。 同 时 
可 以 发 现 ,尽管 多 参数 识别 问题 解 不 唯一 ,但 可 以 通 
过 参数 识别 算法 得 到 参数 的 近似 取 值 ,如 a 近似 于 
0.6, 1A(puTh) 近似 于 0.04 等 ,这 对 于 数学 模型 的 
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实际 应 用 亦 具 有 重要 价值 。 
3.2 FPT 模型 拟 合 及 参数 识别 


利用 FPT 模型 拟 合 混凝土 的 徐 变 实验 数据 ,并 
确定 模型 参数 。 在 贝 叶 斯 算法 中 ,参数 初始 值 % = 
(aa 272i) 取 为 (0.5,1,1,1) ,在 布谷 鸟 算法 
中 ,参数 初 值 在 a2 e [0,1] ,mi e [0,25] ,Hz e 
[0,5] ,ri e [0,100] 范围 内 随机 生成 ,算法 参数 B 
取 为 1.5。 

图 7 为 利用 贝 叶 斯 算法 所 得 的 FPT 模型 中 参数 
oa 、1A(jwaiT%  ) 、1M 的 迭代 结果 图 。 图 8 为 利用 
布谷 鸟 搜索 算法 所 得 的 参数 迭代 结果 图 。 由 图 可 
见 ,在 迭代 一 定 次 数 后 ,两 种 算法 所 得 的 参数 值 亦 趋 
于 稳定 ,说 明 在 FPT 模型 的 参数 识别 问题 中 , 贝 叶 斯 
算法 和 布谷 乌 搜 索 算法 都 是 收敛 的 。 同 样 ,布谷 乌 
搜索 算法 达到 稳定 值 所 需 的 迭代 次 数 更 少 ,说 明 布 
谷 鸟 搜索 算法 在 该 分 数 阶 参数 识别 问题 中 效率 更 
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Fig.7 The Bayesian algorithm iteration results of FPT model 
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8 FFT 模型 布谷 岛 算 法 迭代 结果 
Fig.8 The cuckoo search algorithm iteration results of FPT model 
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图 9 为 FPT 模型 中 布谷 乌 算 法 的 适应 度 迭 代 
图 。 由 图 可 见 , 随 着 迭代 次 数 增加 ,误差 逐渐 减 小 且 
趋 于 一 个 稳定 值 ,进一步 说 明 布 谷 鸟 搜索 算法 在 分 
数 阶 模型 的 多 参数 识别 问题 中 的 有 效 性 。 
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9 ”FPT 模型 布谷 鸟 算法 适应 度 迭 代 图 


Fig.9 The fitness iteration results of cuckoo 


search algorithm for FPT model 

图 10 表示 分 别 将 两 种 算法 所 得 参数 带 入 FPT 
模型 得 到 的 实验 数据 拟 合 图 。 从 拟 合 曲线 可 以 看 
出 ,两 条 曲线 与 实验 数据 测量 值 几乎 重合 ,说 明 贝 叶 
斯 算法 与 布谷 乌 算 法 所 得 的 参数 估计 结果 均 是 有 效 
的 ,在 利用 FPT 模型 拟 合 混凝土 的 徐 变 实验 数据 时 ， 
相应 的 多 参数 识别 问题 亦 存在 解 不 唯一 的 情形 。 在 
实际 问题 中 ,需要 根据 参数 的 实际 意义 与 充分 的 测 
量 数据 ,进一步 确定 模型 参数 的 精确 取 值 。 
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Fig. 10 Experimental data fitting of FPT model 
表 2 为 FPT 模型 的 参数 识别 结果 。 由 表 2 可 
知 ,虽然 两 种 算法 所 得 的 参数 结果 存在 偏差 ,但 与 实 
验 数据 的 拟 合 误差 均 较 小 ,考虑 到 布谷 鸟 搜 索 算法 
的 效率 更 高 ,速度 更 快 , 且 误 差 更 小 ,说 明 相 较 于 贝 


第 40 卷 
叶 斯 算法 ,布谷 马 搜索 算法 在 FPT 模型 的 参数 识别 
问题 中 效果 更 好 。 


表 2 FPT 模型 参数 识别 结果 
Tab.2 The parameter estimation results of FPT model 


参数 贝 叶 斯 算法 布谷 鸟 算法 
a 0.55 0.60 

1/ (p17 ) 0.0375 0.0955 
1/1 0.068 0.0477 

e 0.000 423 0. 000 374 


此 外 ,对 比 表 1、 表 2 可 得 ,用 修正 的 分 数 阶 
Maxwell 模型 与 FPT 模型 所 得 的 拟 合 误差 均 达 到 
10” 量 级 , 且 误 差 都 在 0.0004 左右 ,说 明 本 研究 所 
提出 的 修正 的 分 数 阶 Maxwell 模型 与 FPT 模型 在 模 
拟 混凝土 的 徐 变 特性 方面 均 是 有 效 可 行 的 。 


4 结 论 


本 研究 对 分 数 阶 模型 模拟 混凝土 徐 变 特性 的 多 
参数 识别 问题 进行 了 研究 。 在 模型 方面 ,提出 了 利 
用 修正 的 分 数 阶 Maxwell 模型 与 分 数 阶 Poynting- 
Thomson 模型 模拟 混凝土 的 徐 变 特性 ,通过 与 实验 
数据 拟 合 ,说 明了 两 种 模型 在 模拟 混凝土 的 徐 变 特 
性 方面 均 是 有 效 的。 在 多 参数 识别 方面 ,采用 了 贝 
叶 斯 算法 和 布谷 鸟 搜 索 算法 两 种 算法 同时 确定 了 模 
型 中 的 多 个 参数 ,通过 参数 迭代 图 和 数据 拟 合 误差 
的 对 比分 析 , 验 证 了 贝 叶 斯 算法 和 布谷 乌 搜 索 算 法 
所 得 的 参数 估计 结果 均 是 有 效 的 ,但 相 较 于 贝 叶 斯 
算法 ,布谷 鸟 搜索 算法 搜索 速度 更 快 , 效 率 更 高 误 
差 更 小 ,在 分 数 阶 模型 的 多 参数 识别 问题 中 性 能 
更 优 。 
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